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0\ ; Abstract 

The effective field theory (EFT) for a toy model of nucleons interacting via 
a short range and a long range Yukawa potential is presented. The scattering 
CN \ amplitude in the ^Sq channel is calculated up to next-to-next-to-leading or- 

der (NNLO) using the Kaplan-Savage- Wise (KSW) power counting scheme. 
C^ ' A particular expansion of the amplitude about the pole at low (imaginary) 

J^ ■ momentum is used to derive matching conditions for the EFT couplings. In 

t~^ , addition, after imposing constraints from a Renormalization Group flow anal- 

^^ \ ysis, there are no new free parameter in the amplitude at NNLO. Comparing 

f^ ' the NNLO phase shift to the full model we show that the EFT expansion is 

0^ 



X 



converging. 



^ I. INTRODUCTION 

rH \ Effective field theory (EFT) can be a very useful tool in the study of nuclear physics at 

low energies [|l|- |2^. Most importantly, it provides a perturbative description of physical 



processes at momentum p well below a certain cutoff A. If fields for the lighter degrees 
of freedom are included explicitly, and the effects of the heavier particles are encoded in a 
5^ \ derivative expansion of local operators, one expects effects of the high-energy physics to be 

suppressed by powers of p/A. Furthermore, a perturbative expansion allows a systematic 
estimation of errors at any order of the perturbation. 

The contact terms in the Lagrangian that replace high-energy, short-range interactions, 
have dimensionful couplings. The mass scales that enter these couplings originate in the 
microscopic theory, and they are usually assumed to be "natural", i.e. of C(A). For natural 
sized couplings, the derivative expansion of contact terms in the Lagrangian is an expansion 
in powers of p/A. 

However, constructing an EFT for nucleons requires special care. The reason is the 
appearance of other length scales, in particular the ^5*0 scattering length, a, which is much 
larger than the range of the interactions (a ~ 1/8 MeV ^ 1/m^, 1/mp). In this case the 
coupling constants need not follow the natural scaling of 1/A, but rather may make up the 
dimensions by factors of a. 

To understand the problem better, one first considers nucleons interacting through a 
single short range potential P] . Lacking other "light" scales in the problem, one expects the 



scattering amplitude at low energies to be a simple Taylor expansion in p/A of the exact 
amplitude. The scattering amplitude is related to the phase shift in the ^Sq channel through 

^ = - — \ — , (1.1) 

M p cot — ip 



where p = yME is the momentum of each nucleon in the center of mass coordinates and 5 
is the ^5*0 phase shift. The problem can be analyzed with the tools of quantum mechanics 
(which applies to the non-relativistic, heavy nucleons) where one writes what is known as 
the Effective Range Expansion (ERE): 

1 1 °° 
pcot5 = — + -^r„p2(n+i)_ (^2) 

'^ ^ n=0 

The scattering length a can be arbitrarily large, while the other length scales tq, ri, . . . are 
assumed to all have natural sizes, r„ ~ 1/A*^^""'"^^ p is taken to be small compared to A, but 
not compared to 1/a <^ A. 

Substituting Eq. (|1.2|) in Eq. ( |1.1|) , it can be seen that there is a pole in the amplitude 
at p ~ i/a, which sets a small radius of convergence for the naive Taylor expansion about 
p = 0. This explains the breakdown of the simple approach at momenta p ~ 1/a, and also 
suggests the solution - the EFT must be constructed in such a way as to include the pole in 
the amplitude at all orders. 

In the following expansion of the amplitude was presented: 






M (l/a + ip) 



ro/2 2 , (^o/2) 
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1 + 771 ^TP + 771 ~T^P 

[l/a + ip) (l/a + tp)'^ 
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{l/a + ipY (l/a + ip) 



(1.3) 

In this paper we examine a different expansion I^Gj in which the exact pole in the amplitude 
appears explicitly. We rewrite (|1.2|) as 



1 oo , 

pcot6 = ~-f + -Y.Sn{p^ + l^y . (1.4) 

^ n=0 

where ij is the pole in the amplitude. We expand the amplitude in p/A and 7/A. The coef- 
ficients Sn and r„ are the same at leading order (LO). Comparing the expansions Eqns. ( |1.2|) 
and ( |1.4|) , one can solve perturbatively for the relation between the coefficients r„ and Sn- 
The first two terms are 

1=1- ll'so + 0(7') 

ro = 50 + 272^1 + 0(7'). (1.5) 

The difference between 7 and 1/a enters at next-to-leading order (NLO), whereas the dif- 
ference between sq and tq is only important at orders higher than next-to-next-to-leading 



order (NNLO). Combining Eqns. ( p..4| ) and (|1 . 1| ) gives 
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The EFT for nucleons with exphcit pions has aheady been used P to calculate the ^5*0 
phase shift at LO and NLO. An NNLO calculation is needed primarily to check the conver- 
gence of the EFT expansion. More generally it serves as a probe for better understanding 
of different aspects of the theory, like the appropriate fitting procedure and the role of the 
renormalization group flow. 

The scattering amplitude at NNLO gets contributions from instantaneous pion exchanges 
( "potential pions" ), as in NLO, and also from retardation effects in pion exchange ( "radiation 
pions") and relativistic corrections. The relativistic corrections are very small compared to 
the potential pion contribution even at high momentum. The radiation pion effects have 
been found to be highly suppressed due to Wigner's SU{A) spin-isospin symmetry in the 
large scattering length limit in the singlet and triplet channel of N-N scattering | |27| ]. 

In this paper we analyze a toy model, with some realistic features, but without retarda- 
tion pion effects and relativistic corrections. The EFT for nucleons using the power counting 
scheme introduced by Kaplan, Savage and Wise (KSW), is not completely understood yet, 
and it is beneficial to probe it further with an NNLO calculation in a "cleaner" setup. It 
should be noted that the EFT amplitude for this toy model and the one appropriate for the 
real ^5*0 case differ only by the small relativistic and radiation pions effects. The contribu- 
tion of potential pions to the real amplitude is identical to the amplitude calculated below, 
in Sections ^ and ^U . 



We calculate the ^Sq phase shift up to NNLO in a two- Yukawa toy model which is 
presented in Section |l|, with the corresponding effective Lagrangian described in Section III. 
A brief review of the KSW power counting is presented in Section |V[ Section includes the 
calculation in the effective theory without pions, whereas the NLO and NNLO scattering 
amplitudes in the theory with pions are calculated in Sections [V^ and |V11| respectively. In 
Subsection [VIII A| we derive constraints on some of the EFT couplings. The renormalization 
group flow is discussed in Subsection |V111B| , where the rest of the couplings at NNLO are 
fixed. The comparison with the two- Yukawa phase shift is presented in Section p(| , followed 
by summary and conclusions in Section ^. 



II. THE TOY MODEL 



In this section a theory which can be solved exactly is described. Further below we model 
the low energy behavior of this theory by an EFT. 



We consider non-relativistic nucleons interacting via two instantaneous Yukawa poten- 
tials: 

(r is the distance between the nucleons). This model has several appealing features. It is 
complex enough in the sense that the interactions involve two well separated scales as in the 
real world. In addition, the Yukawa potentials describe the instantaneous exchange of the 
71 and p mesons (without the complications of retardation effects p7| ). 

The realistic ^5*0 one pion exchange amplitude (here given in momentum space) is 

( 9a\ p" ^ 9a 9a ^l .2 2^ 

1^ 2/2^^2+^2 2/2 -^2/^2 + ^2' I • ) 

where the first term is a contact interaction, and the second is the Yukawa part. The contact 
interaction is not included in the toy model, g-,,: is chosen to reproduce the Yukawa part. 
Comparing Eq. (|2.2| ) and Eq. ( |2.1| ), we set 

9. = ^y^. (2.3) 

We take the nucleon mass to be M = 940 MeV, as well as itlt^ = 140 MeV, nip = 
770 MeV, fvr = 132 MeV and qa = 1.25. Qp is tuned to give a large scattering length 
[qp = 13.5 — i> a = —24.7 fm) by numerically solving the Schroedinger equation with this 
potential. The toy model is used in generating all the ^5*0 phase shifts. 



III. THE EFFECTIVE LAGRANGIAN 

The effective Lagrangian for the toy model is described in this section. It is appropriate 
for center of mass energies much smaller than nip. The effects of the short range force are 
reproduced in the EFT by contact interactions. As noted in [^, it is convenient to separate 
the effective Lagrangian in terms of the number of nucleon fields present: 

C = Ci + C2 + -, 

where £„ describes the propagation or scattering of n nucleons. £i is the usual kinetic 
energy term [p|,p8| 

£i = N^ido + W^/2M)N, (3.1) 

where the isodoublet field N represents the nucleons: 

N = 

Next, we write £2 ° ; the part of £2 that has non- vanishing matrix elements between ^Sq 
states. The matrices Pk [@ are used to project onto the ^5*0 state. 




Pk = ^0^2 ®r2rfc, 

where the a matrices act on the nucleon spin space and the r matrices act on the nucleon 
isospin space. The result is 0,H: 

4''^ = -(f )'~''(^o + D^ml + D,ml + ■ ■ ■){N^ P,N)\N^ P,N) 

^%"\\C, + i/^f m^ + . . .)[{N^m)\N^n{vrN) + h.c] 

+ ■■■ (3.2) 

(summation over the repeated isospin index i is imphed.) 

We use dimensional regularization where the Lagrangian is given in D space-time dimen- 
sions and /i/2 is an arbitrary mass scale introduced so that the couplings C2n, D2n have the 
same units in any dimension D. C2n are the coefficients of contact interactions containing 
2n derivatives and D2n ^^^ coefficients of contact interactions involving 2m derivatives and 
2{n — m) powers of the pion mass. We use the convention D2n = D2n-i ^in = D2n ■ The 
ellipses indicate operators with higher powers of derivatives or pion mass insertions. 

For n > (m > 0), the couplings C2n{D2n) multiply a linear combination of operators 
. These operators differ in the way the derivatives act on the legs entering the vertex. In 
many cases factors of internal loop momentum q result in factors of external momentum », 
and there is no need to distinguish between the different parts that make up a C2n{D2n) 
vertex. At NNLO the ffist exception to this rule is encountered. We come back to this point 
in Section IVllB. 



The interaction with the pions is described by a non-local term in the action: 

n J. „— m^|x— y| 

S^ = j d'^x d^y (N^{x)P.N{x)yg f _ 5{x' - /) (n^ {y)P.N{y)) 



47r|x — y| 
d^p dy (N^{p)PMp)y ^2] tI 2 K(p')^.Ar(p')) . (3.3) 



Once the effective Lagrangian has been described, one needs a power counting scheme, 
i.e. a set of rules that determine the relative sizes of diagrams contributing to the scattering 
process. 

IV. THE KAPLAN-SAVAGE- WISE POWER COUNTING 

The KSW power counting scheme was introduced in |^, and is designed to incorporate 
the effects of an unnaturally large scattering length. 

The two versions of the ERE, Eq. ( |1.2|) and Eq. ( [1.4|) , lead to two different expansions of 
the amplitude: Eq. ( |1.3| ) and Eq. ( |1.6D . Both expansions coincide in the infinite scattering 
length limit, but they are not identical. In the next section, the effective theory without pions 
is revisited, mostly to demonstrate some of the consequences of using expansion Eq. ( [I.6[ ). 
In this section, however, we concentrate on the KSW power counting scheme which follows 
from either expansion. 



i^.-i = 




It has been shown [jT| that the leading term in expansion Eq. (1^) is obtained by summing 
up the bubble chains in Fig. |I] with insertions of Cq, giving the LO amplitude: 
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where L is the loop integral 
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(4.2) 



A subtraction scheme that reproduces the KSW power counting is the power divergence 
subtraction (PDS) scheme that was introduced in 0. In PDS one subtracts the poles of L 
in both D = 4 and D = 3 dimension (which gives Xsub = // in Eq. (|4.2| )). Setting /i ~ p, 
all the diagrams contributing to the bubble chains are rendered equal in size (~ 1/p). This 
justifies summing up the bubble graphs at LO. 

In the theory with pions, pion effects are assumed to be perturbative. From the PDS 
scheme follows a power counting scheme - KSW - in which: 

• The expansion parameter is Q/A, where p,^,mT^ ~ Q. 

• A = EZ-l An, An ~ g". 

• The renormalized coupling Cq = Co(/i) scales as l//i ~ 1/Q, and more generally the 
renormalized couplings C2n, -D2J" scale as l/Q"^^- 

Some immediate consequences are: 

• The loop integral L scales like p or /i, i.e. 0{Q). Therefore, adding a Cq and a loop L 
to a diagram does not change its order since the powers of Q cancel. It is now required 
to dress each diagram in the theory by Cq to all orders. (This is diagrammatically 
represented by the "blob", Fig. ^. 

• When the nucleons are put on shell, nucleon propagators scale as 1/Q^ and loop 
integrals J dq'^ scale as Q^. 






Cq Cq Cq 




FIG. 2. Definition of the "blob" 

This theory should correctly describe the physics at momentuin below A ~ rnp/2, which 
is the scale associated with the cut in the amplitude due to p exchange. A calculation to 
order Q" is expected to be accurate up to an error of ~ (Q/A)""*"^. 

In the last three sections, all the components of the computational framework of the EFT 
for nucleons have been identified: the effective Lagrangian (Feynman rules), the regulariza- 
tion and subtraction methods, and the power counting scheme that organizes the diagrams 
of the theory in a perturbative series. In the following sections we apply this formalism to 
the calculation of the ^5*0 scattering amplitude. 

V. THE EFFECTIVE THEORY WITHOUT PIONS 

At energies low compared to the pion mass, even the long range pion exchange is ef- 
fectively point-like. Both the real N-N scattering process and the one in the two- Yukawa 
theory are modeled by the same EFT without pions at these energies, except for relativistic 
corrections. The theory without pions has been solved before 0, using expansion Eq. ( p.. 31) . 
However, since we use the Q expansion Eq. ( |1.6|) , it is instructive to consider the effective 
theory without pions again. 

Expansion Eq. (L3) is not well behaved near the pole, at momentum p ~ i/a. This is 



because higher orders in the expansion contain terms more singular than the leading order 
term. This means that the Q expansion breaks down close to the pole. Expansion Eq. ( |1.6|) 
avoids this problem. All terms have only a simple pole and the expansion is not expected 
to fail even at very low p. 

The structure of the effective Lagrangian when pion effects are integrated out can be 
obtained from the Lagrangian that includes pions by dropping all the D couplings, as well 
as the interaction part 3,^. Note that this only gives the structure of the Lagrangian, but 
the meaning of the C couplings is different in the two EFTs, and the relation between them 
can only be seen by comparing the theories in a region where they are both valid. 

Expansion Eq. ( |1.6D is the exact description of the phase shift at low energies. The EFT 
expansion, A = A-i -|- ^o + -^i H — •, should recover it. We equate the two expansions order 
by order and solve for the couplings C2n in terms of the parameters 7, sq, si, . . . Matching 
the LO amplitude with expansion Eq. (p..6|), we find 
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(5.1) 



At NLO, dressing the C2p^ operators with Co bubbles yields 



^0 = -C2P' 



-A-i 



p" 



Co J ii + ipy 

This does not have the right form to reproduce the Q expansion (|1.6|) , which predicts 



Ao 



p2 _^^2 

(7 + ipy 



To match to the Q expansion at NLO, a contact operator similar to Cq is needed. This 
suggests that a sub-leading piece from Cq is required in order to reproduce the low energy 
expansion. More generally, we expand the parameters in powers of Q ||29||: 



Co 
C2 

C4 



Co,-l + Co,o + Co,! + • • • 
^2,-2 + L'2,-1 + ■ ■ ■ 
Ca-3 + C4 2 + ■ ■ ■ 



(5.2) 



where the second subscript denotes the scaling with powers of Q. The Feynman diagrams 
in this theory can be summed to all orders, and the EFT expansion can reproduce exactly 
the ERE [ig. Working up to NNLO we find in the PDS scheme 
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(5.3) 



Expansions Eq. (L6) and Eq. (O) coincide in the limit a — >■ 0x0(7 ^ 0). This is also 
demonstrated in the vanishing of the sub-leading couplings in this limit. The couplings are 
given in terms of the "experimental" quantities 7, So, Si, etc.. At LO Co,_i is determined 
from 7, at NLO Co,o and C2-2 are given in terms of 7 and So, and at NNLO Co,i, C2,_i, 
and C4,_3 are again fixed by 7 and so- From Eq. ( |5.3| ) it can also be seen that only for 
the leading C2nS, the scaling with powers of \i is the same as the scaling with powers of Q. 
This demonstrates explicitly how knowledge of the running couplings may not be enough 
to determine the power counting of an operator. Here, it is the expansion about the exact 
pole that has lead to this behavior. 

To relate the scattering amplitude to phase shift, we use the relation 
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FIG. 3. The Effective Theory Without Pions. (a) 6{degrees) vs. p(MeV) (b) log-log plot of 
A = I cot Seft — cot 6exact\ VS. p(MeV). The exact phase shift is described by the solid curve. The 
long-dashed curves denote the LO phase shift and error, the dashed curves are the NLO results, 
and the dot-dashed curves describe the NNLO results. 
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The comparison with the exact phase shift of the two- Yukawa model up to NNLO is 
shown in Fig. |] (a). The theory describes the "data" well up to p ~ 50 MeV. To investigate 
the cutoff A of this theory we look at the error plots [^ Fig. |^ (6). We see that at 
p ~ 200 MeV, the errors from different orders in Q expansion become comparable. The 
conditions Eq. ( |5.3| ) guarantee that the expansion Eq. ( |1.6| ) is recovered. Since there are no 
contributions to pcoiS at odd powers of Q in this expansion, the NLO and NNLO curves 
in Fig. ^(6) are the same. 

Fig. D is produced using only the ERE. The EFT approach provides a broad formalism 
in which the pions are naturally included. The next two sections describe the calculation of 
the scattering amplitude to NLO and NNLO in the effective theory with pions. 
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FIG. 4. A 
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VI. THE EFFECTIVE THEORY WITH PIONS: NLO 



In the theory with truly dynamical pions, the contact piece of the pion exchange can 
be reabsorbed in the definition of Co,o- Hence, at NLO, there is no difference in the ampli- 
tude for this toy model and the real N-N scattering amplitude as the radiation pion effects 
and relativistic corrections appear at NNLO. In p[, the ^Sq phase shift has already been 
calculated to NLO. We include the calculation here for clarity. 

At NLO, we need to consider local operators that make up the following four-nucleon 
vertices, classified here according to their Q counting: 



^2,-2 



Only the linear combination Cq^ + m^D2,-2 appears in the amplitude, and it is denoted by 
-Bq.o- From this point on we adopt the convention that the second subscript is omitted from 
leading couplings, so that, for instance, Cq stands for Co,_i, C2 is really C2-2, etc. Also 
note that at NLO, the operator €2^^ can be replaced by C2p^. As explained in Appendix 
B, this is not always the case. 

The diagrams contributing to the NLO amplitude are shown in Fig. ^ The amplitude is 



(Ott) 
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where ^n '" 



Ao = Ar + A}; 
includes all the diagrams with no pion exchange, and Aq 



contains all the 

diagrams with a single pion exchange. It is convenient to express Fig. ^b) and Fig. |^(c) in 
terms of the loop integrals Pi and P2, so that 
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Eq. ( p. 3D , where only the poles have been subtracted, differs from an earlier calculation 
of P2 [0 that involved finite subtractions. The calculation of P2 in coordinate space is 
presented in Appendix A. 



VII. THE EFFECTIVE THEORY WITH PIONS: NNLO 

At this order, there are insertions of three kinds of local operators, listed here according 
to their Q counting: 

0{Q-^) : -Co 

O(Q0) : _5o,o ; VC2 

0(Qi) : -i?o,i ; Vi?2,-i ; -p'C4 

There are 3 new couplings that enter at this order: i?o,i = Co,i + m'lD2-i + mlD^, -62,-1 = 
C2,-i + mlDf^ and C4. 
The NNLO amplitude is 
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FIG. 5. Al 
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FIG. 6. ^i 
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A. A'; 

The diagrams in Fig. ^ are generated by two insertions of 0{Q^) operators or a single 
insertion of 0{Q^) operator, and give: 
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(1 + U_iL)'(-C4P^ - B2,-1P' - Bo,l] 

+z{l + iA-iLfL{-C2p' - Bo,of. 



(7.1) 



B. ^S") 

All the diagrams of Fig. |^ are generated by insertions of local operators of order 0{Q^) 
and a single pion exchange. 

It is here that we encounter for the first time the need to distinguish between the op- 
erators that make up the contact vertices. A detailed discussion is presented in Appendix 
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B. The result is the introduction of two new integrals, Qi and Q2- The one pion exchange 
contribution at this order is 



2(1 + M_iL) ( -^ ) (Qi + QaM-i)), (7.2) 



where 
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and Pi, 2 are defined in Section [V^. 

c. ^5"") 

Fig. 1^ contains the diagrams of this set. The graphs Fig. |^(e), (d) and (f ) can be written 
in terms of analytic functions Ji (p) , I2 {p) and I3 {p) respectively which are presented in some 
detail in Appendix C The rest of the diagrams can be written in terms of expressions that 
have already been defined. We find 

4--) = -il, + 2(^_i)/2 + tiA^iYh 

+PM-i)Pi + 2i{A.ifPiP2 - {A.i)'\P2f. (7.5) 

In the last two sections, we have calculated an expression for the LO+NLO+NNLO 
amplitude in terms of six independent free parameters Cq, Bqq, C2, Bqi, -82,-1 and C4. 
When the number of parameters is doubled (from NLO to NNLO) it is hard to consider an 
improvement of the fit as evidence for the convergence of the EFT expansion. Therefore, it is 
desirable to find theoretically viable ways to reduce the number of independent parameters. 
This is the subject of the next section. 

VIII. ANALYSIS 

A. Matching conditions 

In this subsection we introduce a method for deriving conditions that reduce the number 
of free parameters. The guiding principle is simple: the leading couplings are defined in 
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iJ^f = 



-A 



V 





+ 2 




(c) 



+ 2 



Ai 



(d) 



(e) 




FIG. 7. ^i 



(tttt) 



(f) 



such a way that certain features of the low energy phase shift (hke the effective range 
expansion parameters) are reproduced exactly. Conditions are then imposed on the sub- 
leading couplings so that these features are not modified at higher orders. Similar matching 
conditions were introduced in i2Df . 



From Eq. (|1.1|), the Q expansion of pcot 5 in terms of the amplitude A = A-i + ^o + 

Ai-\ , is 



p cot 6 = ip + 



Att 



An Ao An 

+ 



A^ 



A^ 



MA.I MA-iA-i MA-i\A\ A^i^ 
At low energies, the same quantity is given by expansion ( |1.4| ): 

pcot6 



(8.1) 



-7 + 2^0 (p' + 7') + 2^1 {P^ + 7')' + ■ 



(8.2) 



where pcot 6 has been expanded about p"^ = —7^ (or p = ±27). 

The LO amplitude contributes only a constant to pcot6. Cq is determined by requiring 



ip + 



Att 



MA^i 



-7. 



(8.3) 



This puts the pole in the amplitude at p = i'j. At NLO and NNLO there are more constant 
contributions to p cot 6. The requirements that the position of the pole does not shift at 
NLO, nor at NNLO, reduce to: 
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NLO: 


Ao\p=-i'y — 


NNLO: 


•^1 \p=—i-y U. 



(8.4) 

These conditions determine the couphngs i?o,o and i?o,i- Similarly, one requires that Sq in 
Eq. (^ is determined at NLO by C2: 

df [-ma7,a7J "^^--^ = r°' ^'-'^ 

and that higher order (NNLO) effects do not change sq: 

d 

NNLO: —Ai\„=-i^ = 0. (8.6) 

op 



This condition determines -82,-1 • (Eqns. (|8.4|) and ( ^.61) can be more directly obtained by 
considering the expansion of the amplitude Eq. (|1.6| ).) 

Note that at NNLO, fixing 7 is not exactly the same as fitting the scattering length a 
whereas reproducing sq is essentially identical to fixing the effective range, tq. 

Implementing the matching conditions Eqns ( ^.4|) and ( ^.6|) , introduces a non-analytic 
ttIt, dependence to -Bq.O) -So,i and -82,-1, e.g. 






Further discussion of this point is presented in the next subsection. 

The values of 7 and sq are used to determine Co and C2. One can determine C4 from si, 
the next free parameter in the expansion Eq. ( |8.2| ). However, as in the theory without pions 
(Eq. ( ^.3|) ), it can be shown, using RG analysis, that at NNLO, C4 is completely determined. 
This is done below. 



B. Renormalization Group Flow 

In this subsection, RG flow equations are analyzed. We discuss how constants with 
Q scaling, namely 7 and ?7i^, enter the couplings when the theory with pions is made to 
reproduce expansion Eq. ( |1.4| ). 

The RG equations are derived by requiring that the amplitude be exactly fi independent 
at each order in Q for any value of p. The solutions of the RG equations for Cq and C2 are: 

C,(,).fe(^)\ (8.8) 

The integration constants ^Co and C,C2 are determined by boundary conditions. The other 
coupling constants, denoted here by X, have the general solution 
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X(^^)=^^(M^d\ +/^(Co(/.),i?o,o(/x),C2(^);/i) (8.9) 

where fx is known, and the only freedom in X is in the constant C,x, determined by boundary 
conditions. Specifically, 



Examination of the solutions in Eq. (^.10|) shows that the /i scaling of some couplings 



does not agree with their expected Q counting. To determine the Q counting one also needs 
to know the Q counting of the RG constants C,x- One way of examining this without any 
further assumptions is by finding what values of all six parameters give reasonable fits to 
the data. Following this procedure we have found that at least some of the integration 
constants ^x must have non-trivial Q counting in order to reproduce the low energy phase 
shift (note that for determining the Q counting, the factors of m,^ in the definitions of the B 
couplings must be taken into account). This result is independent of the specific form of the 
low energy expansion used or the particular matching conditions imposed. Using a specific 
fitting procedure, we have found explicit examples of 7 and ?7i,r entering the couplings. 
That 7 can make up Q counting was already shown in the theory without pions, where, for 
example, Cq^ ~ 7^/(Ai — 7)^ (see Eqns. (p.3|)). 

A remark about the m^ dependence of the theory is appropriate. In the usual approach 
to EFT, the effects encapsulated in the contact interactions are the result of integrating out 
fields associated with heavy particles. In the case of the two Yukawa toy model, this would 
mean that all the couplings - Cs and Ds - depend only on the details of the physics at the 
scale rrip where the pion is practically massless, and one would therefore expect them to be 
independent of m,r- More specifically, one expects Cq to be a function only of the pole in 
the amplitude of the high energy theory, that does not include the pions, Cq = Co(7p). 

It is a key feature of the KSW power counting scheme, however, that at LO the amplitude 
has the right pole, or the correct scattering length. This implies Cq-i = Co,-i(7). Note 
that in general the exact pole in the amplitude has a complicated dependence on pion 
physics, and in particular 7 = 7(771^). At higher orders there are contributions from pion 
graphs that move the position of the pole in the amplitude, but these effects have already 
been accounted for. Counterterms must be introduced to cancel these effects, which are 
in general complicated (non- analytic) functions of m-Tr. This has already been shown to 
happen: satisfying the matching conditions Eq. ( |8^ ) (introduced indeed to keep the pole 
from moving), implied explicit m^r dependence of Bqq (Eq. ( ^.7|) ). 

Once it is understood that RG constants such as m^ and 7 can enter the couplings, the 
fi dependence is not enough for determining the Q counting. It is possible, however, that a 
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FIG. 8. (a) The phase shift 6 for the ^^o channel, (b) Log-log plot of A = | cot 6 eft — cot 6exact\ 
vs. p (MeV). The exact phase shift is described by the solid curve. The long-dashed curves denote 
the LO phase shift and error, the dashed curves are the NLO results, and the dot-dashed curves 
describe NNLO results. 

more careful statement is true, concerning the leading couplings. We demonstrate this point 
by discussing directly the coupling C4, which has power counting of 1/Q^- The second term 
in the solution of the RG equation for C4 (Eq. ( p.lO|) ) is ~ ^d/f^"^- In order for it to contribute 
at NNLO, C,Ci should be 0{1/Q). Since C4 has a sensible large scattering length (7 -^ 0) 
limit, the possibility ^(74 ~ 1/7 is ruled out. It is also assumed that C4 is fundamentally 
independent of pion physics, and that absorbing pion effects in it (as in going from Co(7p) 
to 6*0-1(7)) introduces explicit rUn dependence only in the sub-leading couplings. These 
considerations allow setting (,04, = at this order, and have C4 be completely determined by 
Co and C2. 

There are in fact no new parameters at NNLO! 

In the next section we check how well the EFT reproduces the phase shift of the toy 
model when the matching conditions derived in the last two subsections are imposed. 

IX. RESULTS 

In this section the EFT expansion is compared to the exact two Yukawa theory. We 
consider how well the "data" is reproduced by the EFT, whether the expansion is converging, 
and what is the breakdown scale. The results are presented in Fig. ^(a). Fig. |^(b) shows the 
normalized errors at each order, A = | cot Seft — cot Sexactl on a log- log scale. Cq is fixed at 
LO by 7, and C2 at NLO by 7 and Sq. By the RG equations, C4 = C|/Co. 

The graphs show a clear improvement of the fit at NLO compared to LO. There is also 
improvement at NNLO, without introducing any new free parameters. The error plot shows 
that the errors at NNLO are smaller than those at NLO. The breakdown scale, where the 
errors become comparable, can be read off of the error plot to be ~ 1000 MeV. This is 
clearly an unrealistic bound since p exchange introduces a cut in the amplitude already at 
mp/2. 

The table below shows the couplings at /i = m.,^. The entries are given in units of fm . 
To make a direct comparison of the couplings easier, we included the momentum factors 
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from the operators C2P^, i?2,-iP^ and C^p'^, evaluated with p 



m^ 



LO 


NLO 


NNLO 


Co 


C2ml 


Bofl 


-Bo,i 


B2-1 ml 


C,ml 


-3.526 


0.06912 


1.412 


0.5620 


0.2886 


-0.001355 



From the table we see a hierarchy in the sizes of the couplings corresponding to the orders 
in the expansion. This is a verification of the power counting scheme. A rough estimate of 
the expansion parameter is Q ~ 771,^/(^^/2) ~ 1/3, which is also suggested by the table. 

X. SUMMARY AND CONCLUSIONS 

In this paper, an NNLO calculation for the ^5*0 channel amplitude of N-N scattering 
for a two Yukawa model was presented. The amplitude is identical to the contribution of 
potential pions to the real ^Sq scattering amplitude. Matching conditions for the amplitude 
at low momentum were derived. These matching conditions and RG constraints were used 
so that no new free parameters are introduced at NNLO. 

We find significant improvement at NNLO over the NLO result. From the hierarchy in 
the sizes of the couplings we conclude that the expansion is converging. The scale where the 
errors from different orders in the expansion become comparable is about 1000 MeV, which 
is higher than the expected breakdown scale for this theory. 

The matching procedure employed in this paper could be applied to the experimental 
N-N scattering data. Since the relativistic corrections and radiation pion effects are small, 
we expect the analysis and observations in this paper to hold even for real nucleons. 
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APPENDICES 

In many cases it is easier to calculate the integrals in this theory in coordinate space. 
Finite integrals of this type are Ji, I2 and J3. When the integrals are divergent in D = 3 
and/or D = 4, some care is needed[|. P2 is such an example. 



^We would like to thank Prof. L.S. Brown for showing us how to extend the integrals off Z) = 4 
dimension, and specifically for providing us with an example the calculation of P2 is based upon. 



The calculation of P2 in coordinate space is presented in Appendix A. In Appendix B 
we discuss the special integrals appearing first at NNLO when a 6*2^^ acts inside a loop 
with a pion across. I3 is calculated in Appendix C, whereas for Ji and I2, only the result is 
presented. 

Appendix A: P2 

After performing the q^ integral and for convenience making the integral dimensionless, 
Eq. ( |6.3| ) gives 



1 1 



X 



P-p^-te (l-q)2 + r 



;io.i) 



where p = p/m.,^ is dimensionless. Fourier transforming to coordinate space yields 

P2 = -^l^^f "' {^^''^^ M^ j d''''^GM{x)GN{-x)G^{-x), (10.2) 

where 



Gn{x) 



(27r)^-i k^-p"^- le 

are the rescaled (p = p/ttLt^) Fourier transformed nucleon and pion propagators in D dimen- 
sions respectively. Now, consider the general form 

/= I d'^-^^WGiK.-x) 

i 

= r(^) (2nf^-'^/' ' ^ U j ^^ ^"'''^ 

= p.D-iN .^ ,3(D-i)/2 y ^^^' ' n('^i^) ' ^iz_3(«ia;), (10.4) 



where 






2 



(27r)^-i F + K 

1 , ^ -D-3 



(27r)^^"^^/^ Va;/ ~ 



— Kd-3 (KiX) 



D^4 e '^'"^ 



Anx 
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10.5) 



and Kij{z) is the modified Bessel function. Tlie divergence in tlie above integral comes from 
small distance x, corresponding to UV divergence. The integral is split into its finite and 
divergent pieces as 

/ = /< + /> 



f7DEi^7^-J3(D^i)72y ^^^^ '^n('«*^) ' ^£_3(/t,a;). (10.6) 



27r(^-i)/2 1 

5 

Using the asymptotic form of the modified Bessel function 

2''-^T(iy) 






we get for / 



< 



J- T-fD-l\ /o.n3(D-1)/2o onii^ ''^ 



r(^) (27r)'(^-^)/' 8 - 2D Y 2 

D^4 1 1 + 1E + ln(47r) + 2 In(^) 

' 327r2(D-4)^ 327r2 

The other integral J^ is finite in D = 4, and gives 



+ 0{{D-4Y), (10.7) 



oo 
<5 


-(K1+K2+K3) 
X3 




7s + ln(Ki + K2 + Ks) + 


HS) 



IGtt^ 
and so, 



;io.8) 



/ = /< + /> = -1 + l-7i. + ln(47r)-21n(., + ., + ^3)_ 

327r2(L)-4) 327r2 ^ ' 

For the integral P2, using 



ni = K2 = y-p^ - ie = -ip; K3 = 1, 
we finally have, after subtracting the pole at Z^ = 4, 

P2 = -^^^'^'^ (-7i^ + 1 + H^) - 2 In(^) - 2 ln(l - 2tp)\ . (10.10) 
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FIG. 9. Diagrams in which one cannot use the replacement V^ 
are (a) Pi and Qi, and (b) P2 and Q2- 



p . The integrals that result 



Appendix B: Qi,Q2 

In this appendix we discuss the action of a 6*2^^ operator in a loop containing pions. 
From the Lagrangian Eq. ( p.2|) we rewrite the operator corresponding to C2 as: 

C2 



'-{Qt+q^: 



(10.11) 



with qi, q2 the single nucleon incoming and outgoing momenta, in the center of mass 
coordinates. At tree level qi = q2 = P, and the contribution from each piece is ^p"^- This 
is, in fact, the contribution from each piece of C2 even when the legs are not external, so 
long as the loops attached to the vertex have no pions. In the diagrams of Fig ^ , however, 
corresponding to the types Fig. H(a) and (d), the situation is different. The vertex reads 
^(p^ + q'^), with q being the loop momentum. Using the fact that the loop still contains a 
nucleon propagator we write (after performing the g° integral): 

.2 1 



Q 



l+p' 



%e 



q^ — p^ — It q'- — p"- 

The second term looks like the p^ term; the contributions (each with C2/2) add up to give 
the "naive" answer (i.e. the answer one gets by replacing V^ -^ p^\ But the "1" remains, 
giving rise to different integrals. 

In the case of Fig. |^(a), when the derivatives act on the external legs, the loop integral 

; when they act on the internal legs, the resulting integral is 



is P\ as defined in Eq. 



2f 



miM ^ 



A-D 



W-l 



q 



(27r)^-i (q - P) 



mt 



PDS 9a 2 



2f 



miM 



/i — m,r 

4:71 



(10.12) 



Similarly the contribution from diagram Fig. |^(b) can be expressed in terms of the loop 
integral P2 as defined in Eq. (|6.3|) and 



Q2 



-t—m^M 
QiL 



PDS 



A 

2p- 



.^lArnlM^i 



/i 
2 

fi- 



8-2D 



q 



:D-1 



1 



1 



1 



(27r)^-i (27r)^-i P - p"^ - it (q - 1)2 + ml 



A-K 



m fi + ip^ 



4:71 



(10.13) 
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In carrying out the subtractions in Qi and (52, some care is needed. The poles and the 
corresponding counterterms should be calculated treating Cq perturbatively. 

Appendix C: Ii,l2,Iz 
From Fig. |^(e), we define 

d^k 1 ,„ .2 



J f27r)3 p - „2 _ ,/(k' P) 



(27r)'^ k"^ — p"^ — ie 
J(k, p) = / ^ , \ . 7 r^— ^- (10.14) 

/(k, p) is finite in D = 4 and one can therefore use the simple form of the propagators in 



coordinate space ,Eq. ( 10.5 ), which gives 



and so, 



/(k, p) = 2^ j dx^-^^^e-^^--'">, (10.15) 





1 /• „ , , e-(™-"*P)(^'+2') sin(xA;) sin(|/fc) 
dkdxdy- 



32tt^ J xy k"^ — p"^ — ie 



1 7,/Lflr!;: L,!!£M£!£(Mi^ ^,,,,^ 



IQti^ J J stil -t) J k^-p'^-ie 

^ ^ 

For p > 0, performing a contour integral in k gives 

■i r 1 ^ „-i2pst _ 1 



QAtt^P J t{l~t)J s 



QATT^prriT, \ \l — i2pj \l — i2pj 

+ Li2{l - ip) - Li2{l - i2p)) (10.17) 



where p = p/m.,^ and the dilogarithm, Li2{z), is defined as 

^ln(l-t) 



Li2{z) = — dt 



/i and I2 are determined from Fig. |^ (e), (d): 



t 
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^1 = M^1 <^^1 



'' ^ I 



d^q 1 



(27r)3 q^ — p^ — ie (q — p)^ + m^ (q — p')2 + m^ ^ 
^ '^L(p) + M(p)) 



327rm^p^ 

'^' 1 , /(l + 4p2t)(l + 4p2(l-t)) 



L(p) = / dt- In 

^'^ J tl-t) 



„ t(l-t) V 1 + V 



1/2 

M{p) =2 f dt- (tan-\2p^Jl + 4:pH{l - t)) - tan-\2p) 

^ ^ 

= ^TT^-iln^ (2) + ^ In (4) In (1 - 2 i p) + 2 TT In (2 - 2 i p) 



In (2 (1 + i p)^ (1 - 2 i p)) In (1 + 2 i p) - tan^^ (2 p) In ^ 



4(l + p2 



+ 4p2 



-2 tan"^(p) ln(l +4p^) -2iLi2 (- + ip 



-2iLi2{-l-2ip) -2iLi2{2-2ip), (10.18) 

where the average is over the outgoing momentum p', and 



2 \ 2 

9 \ ^4n^2] 



^2 = I T^ mlM^h 



I 



2 



2P 
d^q d^l 1 1 1 



(27r)^ {2'kY q^ — p^ — it (q — p)^ + ml P — p^ — ie (q — 1)^ + m. 

i r°° , , (i + x + p\, /(a; + p)2 + l\ 1 

dxm m ' ' 



2 



1287r^m^p j-00 \i — x + pj \{x — pY + 1 j x"^ — p'^ — it 

7r2 + 4i7r ln(2-2?p) + ln(2) ln(l + 2ip) 



128 7r2 p2 m^2 

n4 



+ In (1 — 2 i p) In ( ^^^ —^ — ^—^^ r~T ) ^ -^^2 [7, ~ "^ P 



l-2ipf {l + 2tp)J \2 

+Lt2(- + tp)+4Lt2{2-2tp)Y (10.19) 



The exphcit expressions for Ji, I2 and I3 presented above are vahd only for p > and 
small positive imaginary p (|p| << 1). For other values of p, the integrals have to be done 
taking into account the poles and cuts in the integrand appropriately. 
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